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EXTENDED ABSTRACT 

An empirical comparison is made between three algorithms for computing 
the points of intersection of two planar Bezier curves. The algorithms compared 
are: 

Al. The well known Bezier subdivision algorithm, which is discussed in [Lane’80] 

A2. A subdivision algorithm based on interval analysis due to Koparkar and 
Mudur [Koparkar’83] 

A3. An algorithm due to Sederberg, Anderson and Goldman which reduces the 
problem to one of finding the roots of a univariate polynomial [Seder- 
berg‘84]. 

The details of these three algorithms are presented in their respective references. 
This abstract mentions major implementation choices, and assumes that the 
reader is familiar with the three algorithms. 

The execution times of Al and A2 are highly sensitive to implementation. 
We observed that the decision on when to terminate the subdivision can effect 
execution speed by an order of magnitude. We based our termination criterion 
on a suggestion from [Wang’84] that a curve be subdivided until its deviation 
from a straight line segment is no larger than a value e. If (x,-,y,) i = 0,1 ,..., n are 
the control points of a degree n Bezier curve and 

L o = max (I*,- - 2rr,- +1 + x {+2 \,\ y { - 2y i+1 + y I+2 |), 

0< t< n-2 

then 

v/3n(n-l)L 0 
r ° = log < — & — 

where r 0 is the number of times the curve must be subdivided in order for a 
curve segment to approximate a line segment to within a tolerance e. 
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For Al and A2, the convex hull was taken to be a bounding rectangle. For 
Al, subdivision was performed using the geometric construction algorithm and 
the bounding rectangle was determined by the new Bezier control points. For 
A2, the polynomial equations were converted to standard power basis, and subdi- 
vision was performed by simply evaluating the curve segment at its parametric 
midpoint, using Horner’s algorithm. 

For A3, the implicitization procedure was performed using Bernstein polyno- 
mials. It was noted that standard polynomials introduced excessive roundoff 
error for curves of degree larger than three. A critical part of the implementation 
is a polynomial root finder which finds only real roots. We developed our own 
root finder for this purpose, which usually runs significantly faster than, for 
example, the standard IMSL root finder which finds all roots (real and complex) 
of a polynomial . 

We tested these three algorithms on pairs of Bezier curves of degree three, 
four, and five, with several different numbers of intersection points. Our goal was 
to compute the intersection points to five decimal places of accuracy, and we 
experimented to find the optimal value of r 0 for this purpose. 

Results: The tests were run on a VAXl 1/750, with floating point accelera- 
tor, under UNIX 4.2 BSD. Our test results indicate that for two cubic Bezier 
curves which intersect once, A3 is at least twice as fast as the other two algo- 
rithms. For more than one intersection, A3 does even better — if the two cubic 
curves intersect nine times, A3 is typically five or six times faster than the others. 

For curves of degree four, A2 and A3 are comparable, and Al is the slowest. 
For curves of degree five, A2 is typically twice as fast as the other two algo- 
rithms. 

Multiple intersections were not thoroughly tested. Our implementation of 
Al and A2 did not properly detect the presence of tangency. A3 handled 
tangency with no problem by reporting a multiple polynomial root. 
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